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We present a new method to measure rheological response of liquids confined to nano-scale which 
exhibit a considerable slow-down in dynamics compared to bulk liquids. The method relies on using 
a robust force sensor that has a sensitivity to measure shear forces in the range of 1 to 10 nN and 
a stiffness of 5.5 x 10* N/m to avoid thermal noise. Off-resonance operation ensures application 
of a range of shear frequencies . This range is higher than the inverse of the system's mechanical 
relaxation time and allows the measurement of non-linear effects emerging due to significant "slow 
down" in dynamics upon confinement. The instrument is a step forward towards resolving the 
controversies about the nature of nano-confined liquids. 

I. INTRODUCTION 

The role of nano-confined liquids in many areas ranging from tribology to biology is universally accepted [IHl]- 
Such confinement by solid walls induces significant alteration in their physical properties such as diffusion coefficient 
[6] and viscosity [3 [8], furthermore the liquid possibly undergoes solidification [71 19]. However, the details of this 
change in behaviour have been argued upon ever-since the pioneering experiments that measured oscillatory solvation 
forces on liquids confined to molecular- level thickness [TUl - fT^ . In particular, the measurement of viscosity of confined 
liquid using different experimental approaches have led to contradictory results [131 [T5HT5] . The Surface Force Appa- 
ratus(SFA) measurements have reported that a hydrogen-bonded liquid like water shows no change in viscosity upon 
confinement down to 4 nm and below [17, T5^, whereas other non- associative fiuids such as Octamethylcyclotetrasilox- 
ane(OMCTS) exhibit significant solidification upon confinement [5] . However, many other measurements of shear 
viscosity of nano-confined water using SFA and Atomic Force Microscope(AFM) have reported orders of magnitude 
increase in viscosity under confinement [7} B] • Recently both the shear and normal measurements performed using 
AFM have suggested that the relaxation time of nano-confined water films is orders of magnitude more than that of 
bulk water [H [50] ■ 

The fluid volume involved in the viscosity measurement of nano-confined liquids is of the order of atto-litres to zepto- 
litres, and the shear forces that they generate are of the order of nN. Therefore conventional rheology instruments are 
not capable of revealing the rich non-linear behaviour of these films. SFA and AFM set-ups are modified to measure 
the shear stress developed due to straining of nano-confined thin films. In both SFA and AFM experiments the shear 
frequency is limited from a few Hz to 1 KHz. In order to measure the change in relaxation time upon confinement 
we need a wider shear frequency range, typically larger than the inverse of the system's relaxation time. Also jump- 
to-contact instability of the confining surfaces is a serious issue in controlling the film thickness while shear straining 
these films. It is therefore important to develop new techniques employing force sensors which have large normal 
stiffness to avoid the instability and can perform in wider shear frequency range with force sensitivity of the order 
of nN. This paper reports the development of an instrument for oscillatory rheology on nano-confined liquids with 
unprecedented shear frequencies and near infinite normal stiffness. We also develop the necessary analysis method to 
quantify relaxation time of the nano-confined liquid. 

The new instrument proposed in this paper employs tuning fork as a shear force sensor in rheology of nanoconfined 
liquids. Quartz crystal tuning fork is a commonly available device used as time-keeping element in watches. In recent 
years quartz tuning fork has been used as a robust force sensor in AFM with low signal to noise ratio due to its high 
quality factor of resonance and very high stiffness compared to conventional AFM cantilevers [2TH21]. Tuning fork has 
not been used so far to measure shear response of nano-confined liquids since the charge accumulated due to strain on 
its electrode surface is used to measure the the bending of the arms. This is difficult to measure it liquid environment. 
We use tuning fork as shear stress measuring force sensor by keeping it outside the liquid and mounting a fibre with 
a tip on one of its arms. The tip is immersed in liquid and the tuning fork is used to measure the stress on the fibre 
tip by measuring the difference in current through the arms. The amplitude and phase measurement is related to 
shear stress and the complex damping coefficient by modelling the fibre-tuning fork system as two elastically coupled 
oscillators and solving the resulting differential equation. 
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The paper is organized in the following manner. We first discuss the feasibility of using tuning forks for shear 
measurements in nano-confined water followed by the theory and modelling required to quantify the experimental 
data. We describe the design and development of the instrument and finally present initial results on nano-confined 
water using the instrument. 



II. OFF-RESONANCE TUNING FORK FOR NANO-RHEOLOGY 




FIG. 1: Schematic of the Experimental set-up 



Figure 1 describes the tuning fork as a proposed force sensor in measurement of rheological response of nano- 
confined liquid. A fibre with a sharp tip at its end is attached to the free arm of the tuning fork. It extends further 
from the arm to immerse itself in the liquid. We oscillate one arm of the tuning fork using a piezoelectric crystal. 
The resulting amplitude of the second arm is measured using a lock-in amplifier. The change in this amplitude upon 
bringing the tip close to the substrate placed at the bottom of the liquid cell is a measure of shear stress generated in 
the film; while phase gives us information about the nature of the nano-confined liquid, namely solid-like, liquid-like 
or non-Newtonian fluid-like. 

If the viscosity of the nano-confined liquid remains unchanged compared to its bulk value, the amplitude should not 
change upon approach. If the amplitude reduces upon approach but there is no phase change from the bulk value, then 
the film is more viscous but still Newtonian. A change in both the phase and amplitude suggests a non-Newtonian 
behaviour. Next we discuss the technical aspects of quantifying the experimental data. 

In most applications involving tuning fork as a force sensor, whether AFM imaging HU or force measurements 
during atom manipulations p5] . it is driven on resonance. An atomically sharp tip is mounted on its free arm. It 
is important to keep the tip mass minimal in order to avoid serious loading effects by the tip itself. Secondly, the 
dimensions of the tip are chosen such that the spring constant of the tip is orders of magnitude higher than that of 
the tuning fork and the tip bending is negligible compared to that of the tuning fork arm. The shift in the resonance 
frequency due to interaction between the tip and substrate is measured and is used as feedback control parameter in 
imaging. The shift itself is related to interaction stiffness if the amplitude is small enough compared to the interaction 
length scale ^2 5). However in our application the tip length is significant owing to the fact that the tuning fork remains 
outside the liquid. Secondly, in order to perform shear rheology we need to perform measurements in a range of 
frequency, at least covering one or two decades around the inverse of longest relaxation time the film may possess. 
Therefore we can not restrict ourselves to resonance. In this section we describe the feasibility of using tuning forks 
off resonance with a fibre tip having a lower stiffness than the tuning fork arm itself. 



A. Off-resonance operation 



First we discuss the issue of how to perform measurements in the sub resonance regime. Below resonance the two 
arms of tuning fork move in the same direction and the potential drop across the electrodes is expected to be zero. 
We performed finite element analysis simulation of the tuning fork driven below resonance in the range of 2 to 20 
KHz using COMSOL. Figure 2 shows the bending in the arms when driven off resonance at frequency of 10 KHz. 
Clearly, there is excess bending in the free arm compared to the one mounted on the piezo-crystal and the phase lag 
between the arms is zero. We argue that this extra bending emerges due to inertia. This signal is the measure of 
shear strain amplitude. Experimentally when a tuning fork is driven below resonance, there is indeed a measurable 
signal that arises due to the difference in the current through the two electrodes. Typically with gain of 10^, this 



3 



difference results in lOOmV signal in lock-in amplifier which corresponds to roughly 20nm of free arm amplitude. This 
is discussed in more detail in later sections. 



freq_smpz3d(l)=10000 Subdomain: x-displacement [m] Max: 4.042e-10 




FIG. 2: Finite element analysis Simulation of the tuning fork and piezo arrangement at off-resonance frequency 
using COMSOL.One arm of the fork is oscillated at lOKHz with amplitude of 0.36nm. The amplitude of the free 
arm is more by 10% 



B. Modelling of the shear force sensor 

In oscillatory shear rheology, the phase lag between a sinusoidally varying strain and resulting stress is used for 
measurement of linear viscoelasticity of complex liquids. Typically for purely elastic response the stress is in phase 
with the strain whereas for purely viscous response the phase difference is 90 degrees. For viscoelastic materials the 
phase is between these two extremes . The shear stress is related to shear strain through the relation a = G*7, where 
tJ is shear stress, G* is complex modulus and 7 is shear strain. The storage (elastic) modulus G", and loss (viscous) 
modulus G" is defined as G* = G' + iG" , where i is an imaginary number. The relaxation time of linear viscoelastic 
system is r = {G' /G").l/u! [26]. The complex viscosity is defined a.s ij = rj' — irj" . 

From the foregoing discussion it is clear that in order to measure non-linear behaviour of the nano-confined liquid, 
it is crucial to measure the stress due to applied strain and their phase relationship. The shearing tip experiences a 
stress once it is immersed in bulk liquid. Our goal is to measure the change in this oscillatory stress upon confinement 
of liquid by the tip and its phase relationship with the strain. This is straightforward if we measure the change in 
amplitude of the tip with respect to its position from the surface. The phase can be measured with a lock-in amplifier. 

It is noteworthy that in our measurement scheme the tip extending further from the tuning fork end has a stiffness 
k2 which is much smaller than the tuning fork arm stiffness fci . Ideally, the amplitude and phase at the end of the 
tip should be measured for rheological response of the confined liquid, but experimentally it is not possible. In our 
experiments the amplitude and phase are measured at the tuning fork arm and not at the end of the fibre. In order 
to measure viscoelastic response of nano-confined liquid film from experimentally observed tuning fork amplitude and 
phase, we require the development of a model for the the proposed force sensor. This theoretical treatment is different 
from previous such measurements using AFM and SFA. We model the tuning fork and the attached fibre as a system 
of two coupled springs. Figure 3 describes this modelling in detail. We seek the relationship between the phase at the 
tuning fork arm and the quantities ki, ^2, Li, L2 ,7, where 7 is the damping coefficient. 

The differential equations for the coupled system are 

mix'i —A() sm{ujt) — ki{xi — Li) + k2{x2 — xi — L2) (1) 
1712X2 = - ^X2 - ^2(2:2 ~ xi- L2) (2) 

where Aq is the amplitude of the driving force, ^1,^2 are equilibrium lengths of the two springs. 7 is the damping 
coefficient of the liquid and xi,X2 are the instantaneous tuning fork arm and tip positions respectively. 
Taking the solutions of the form: 

xi = Aie*"* X2 = Aae'"^* 
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FIG. 3: The model of the tuning fork and the fibre. The fibre and the tuning fork arm are treated as springs with 
stiffness ^2 and ki respectively. The damping coefficient 7 is provided by the liquid in which the tip is immersed. Li 
and L2 are the equilibrium lengths of the two springs. 



we get the equations in the matrix form as: 
-miLu'^ + fci + fc2 



-m2ijJ + i^UJ + k2 



j4osin(a;i) + kiLi — k2L2 

k2L2 



Using LU decomposition, we obtain the solutions as: 



Xi 



X2 



Aq sin(wt)(fc2 + *7'^ ^ m2'-o'^) + {i^LO — miuP'){kiLi — fc2-^2) + fcifc2^i 
^2(^1 — miw^) — m2Cj^(fci + fc2 — miw^) + i"fuj{ki + fc2 — miw^) 
A:2j4o sin(a;i) + fcifc2ii + kik2L2 — k2L2miLu^ 
^2(^1 ~ miuj^) — •m2io'^{ki + ^2 — rriiLU^) + ijuj(ki + fc2 — miw^) 



(3) 

(4) 
(5) 



To get the relationship for the phase difference between the oscillator and the tuning fork, as well as that between 
the oscillator and the fibre, we analyse the time dependent part of the solution of xi and X2 ■ 



61 =tan ^ 



- tan 



[(fci + k2 ~ miw2)(72cj2 _ rn2UJ^{k2 - m2Cj2)) + ^2(^1 - miLU^)ik2 - m2UJ^)] 
^Lo{ki + k2 — TOio;^) 



m2UJ^{ki + /c2 ~ TOiCj^) + k2{ki — miio^) 



(6) 
(7) 



We plot 01 and 02 versus damping coefficient 7 in figure 4. The phase of fibre tip, 02 immersed in liquid quickly 
takes the value of 90 degrees for damping coefficient of 0.1. This is a typical value of 7 obtained from water viscosity 
and the area of the fibre that is immersed in the liquid. 0i too converges to degrees. As expected, the tuning fork 
arm which stays out of liquid remains in phase with the driver and the fibre tip immersed in Newtonian viscous liquid 
gives phase difference of 90 degrees. Our solutions for ipi and 02 thus give a familiar result and the solutions we have 
arrived at are consistent with the physical picture of conventional response of Newtonian liquids. 

Solving for 7, equation 6 yields a complex damping coefficient for any value of 0i greater than 0.1 degrees. Figure 
5 shows plots of real and imaginary part of damping coefficient 7 versus 0i . For a fixed geometry of the fibre moving 
inside the liquid, the damping coefficient 7 is proportional to viscosity rj. A complex damping coefficient means a 
complex viscosity. This is a remarkable result, since it allows us to measure the non-Newtonian response of the liquid, 
if any, from the measurement of 0i. Indeed in our experiments when the fibre is lowered into the liquid from air, there 
is no change in 0i. 0i starts to take values in the range of to 5 degrees once the fibre tip reaches within a few nm 
from the substrate, indicating a non-Newtonian response of the nano-confincd liquid. The experimentally measured 
phase difference 0i is used to calculate the real and imaginary part of damping coefficient, namely 7 and 

The complex viscosity of non- Newtonian fluids is defined as 77* — rj —irj and the relaxation time as r = (77 /rj ).l/uj 

m- 

The relaxation time in terms of real and imaginary part of 7 is given by 



T = (7 /7 ) 



1 



(8) 
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FIG. 4: Plots of i/ii and (^2 using equations (6) and (7), the . The fohowing parameter values are used : 
fci = 55090N/m, ki = 140N/ni, mi = 2.2 x lO^^kg, = 2.7 x lO'^kg and uj = 20KHz(a) The phase (jji of the 
tuning fork and (b)the phase of the tip 02 with respect to the driver for different values of 7. Note that for 
Newtonian liquids with reasonable viscosity the 01 converges to zero and the 02 converges to 90 degrees, consistent 
with the ideas in conventional rheology. 




(a) (b) 

FIG. 5: The Real and Imaginary part of damping coefficient 7 with increasing 0i 

III. THE INSTRUMENT 



A. The force sensor 



Tuning forks arc procured from local market having following dimensions for each arm: length L = 3.5mm, width 
w = 350/L(m and thickness t = 700/x m with the electrodes covering the length Lg = 3mm. After carefully removing 
the tuning fork from its protection cap, we mount it on a piezo crystal to drive it mechanically as shown in figure 6. 
On resonance the free arm oscillates with much larger amplitude compared to the one sitting on the piczo-crystal and 
the difference in the charge accumulated due to strain in each arm is measured using a pre-amplifier as shown in the 
circuit diagram. A gain of lO'' is used to convert these oscillations into measurable signal using a lock-in amplifier. 
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FIG. 6: Tuning Fork with fibre tip attached to the free arm 



The resonance characteristic of the tuning fork without loading the fibre tip is measured. A typical quality factor of 
800 is obtained for the operation in ambient environments. A single mode fibre with outer diameter of 127/im and core 
diameter of 9/im is pulled using a CO2 laser based pipette puller to form a sharp tip. The fibre tip is then mounted 
on the free arm of the tuning fork as shown in the figure and the resonance characteristic is measured again to see 
the loading effects. Typically the resonance shift before and after mounting the fibre tip is 1-2 KHz. To estimate the 
amplitude of the free arm in nm from the measured signal in lock-in amplifier we use the geometry of the tuning fork 
and the piezoelectric modulus of quartz following ^21]. The total charge developed per unit deflection of the free arm 
is: 

g/z' = 12d2ifc/ie(^Y-i) A' 

where ^21 is the piezoelectric coupling constant. 
The theoretical sensitivity is given by 

Stheory ~ 27r/ X Rx q/z'. 
For / = 29284iJz, and R = lOAin in our gain circuit 

Stheory = 34.3635 mV/nm. 

B. The coarse and fine approach mechanism 

The instrument is designed and built for shearing the nano-confined liquid between a sharp tip and a flat substrate. 
This requires the actuation of tip towards the substrate with precision of nm and the range of few mm. 




FIG. 7: A photograph showing the X-Y-Z coarse and fine approach mechanism 

For approaching the fibre tip from separations of mm to nm we used the inertial sliding approach mechanism used 
first by D. Pohl ^28^, also discussed in detail in references [211 and [SD]- The entire approach mechanism is placed in a 
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heavy steel base where the central part is carved out to accommodate coarse and fine X-Y-Z positioner. To facilitate 
coarse motion a steel piece called hammer is attached to a piezo tube (hammer piezo). On top of the hammer piezo, a 
quartz glass tube is fitted. This tube is inserted in a connector having a copper leaf spring that holds quartz tube. The 
tension in the leaf spring can be adjusted to serve the controlled vertical sliding of the tube against gravity. The quartz 
tube makes contact with the connector piece along flat walls. Electrical pulses with slowly rising edge and fast drop 
off are applied to the hammer piezo. During the initial slow rising of the pulse, the hammer piezo expands (contracts) 
moving the hammer up (down) and the tube remains clamped at the leaf spring. When the voltage suddenly drops 
to zero, the hammer remains in its place due to inertia pushing the tube up (down) achieving a net movement in 
each pulse. The height and width of the pulse can be controlled to achieve the desired step length, typically tens of 
nm. The connector bearing the quartz tube is fixed onto a coarse X-Y scanner. The quadrant piezo tube (scanner 
piezo) is mounted inside the quartz tube and moves with the hammer piezo, hammer and the quartz tube assembly. 
A liquid cell is magnetically attached on the top of the scanner tube. The cell comprises of two circular Teflon pieces 
clamped together with a Viton 0-ring between them. The top Teflon piece has an opening to allow access to the fibre 
tip mounted on the Tuning Fork. The liquid is filled in the cell from the top opening with a glass syringe. 

The coarse X-Y positioning of the sample also uses inertial sliding principle. The X-Y positioner consists of a heavy 
steel plate called top plate in the figure 7. The top plate has grooves cut out in it for X motion. Similarly, the grooves 
are cut in a perpendicular direction in the instrument base for Y motion. These grooves are seen in Figure 7. Shear 
piezos in alignment with the grooves are glued on both sides of the middle plate. Ruby balls are glued in the recesses 
formed on macor pieces that sit on the shear piezos. For X-Y motion, the slow rising and fast falling pulses are applied 
to shear piezos similar to Z coarse motion described previously. The central part of both the top and middle plates 
has a circular hole for the passage of the Z-coarse positioner and the scanner. The connector holding the Z positioner 
is mounted on the top X-Y plate as shown in Figure 7. All the actuation signals and the feedback loop are generated 
using Lab View. National Instruments data acquisition card (PXT6259) is used for interfacing the instrument with the 
computer. We calibrated the Z-scanner piezo after mounting the liquid cell using a Michelson Morley interferometer. 
This gives the sensitivity of the Z-scanner as 18 nm/V. A circular plate with a hole on the top is used to cover the 
entire approach mechanism assembly. The AFM base is kept on an active vibration isolation table from Accurion, 
micro series model. 



IV. MEASUREMENT METHODOLOGY 

In order to keep the tuning fork out of the liquid we attach a fibre to its free arm so that the fibre tip extends 2 to 
3mm out from the fork and the tip end of the fibre remains immersed in the liquid as shown in the figure 1. Freshly 
cleaved mica is clamped below the Viton O- ring of the liquid cell and the cell is immediately filled with Millipore 
water. The fibre with tip at its end is slowly lowered into the liquid using coarse Z-motion. The change in amplitude 
upon engaging with surface is largest on resonance. We use this fact for engaging the tip with the surface. The tuning 
fork is oscillated on resonance and the tip is approached towards the surface in auto-approach mode. A set of pulses 
are applied to the hammer piezo and between two successive sets the scanner piezo expands its full range to confirm if 
the desired set-point is reached. After reaching the set-point amplitude, the pulses are stopped and the scanner piezo 
is retracted back fully. The whole procedure is repeated with sets having smaller number of pulses and lesser pulse 
height in order to bring the engagement point close to zero Z-voltage on the scanner piezo. This ensures that the tip- 
substrate engagement point is always within reach while acquiring the data. The tuning fork is now oscillated below 
the resonance and the fibre is approached and retracted. The rheological response of the nano-confined liquid between 
the tip and the atomically smooth substrate as the tip is progressively approached towards substrate is inferred by 
measuring the phase lag and amplitude change in the tuning fork arm. Using lock-in amplifier the amplitude and 
phase of the arm of the tuning fork carrying the fibre is measured. The change in amplitude with respect to the 
normal position of the tip from the surface gives us the change in force the tip is experiencing compared to the bulk. 
Figure 8 (a) shows that the amplitude decreases upon approach within about 5 nm from the surface. This reduction 
in amplitude indicates larger drag force upon confinement. Figure 8 (b) shows the phase 4>i versus the separation. 

V. RESULTS AND DISCUSSION 

As discussed in section II, the phase relationship between applied sinusoidal force and the tuning fork's free arm 
can give us the nature of the nano-confined liquid. Any possible phase difference between these two is a indication of 
a non-Newtonian behaviour. Figure 8 shows that there is finite phase change upon approaching the tip within about 
5 nm from the substrate. In AFM experiments, the zero of the tip-substrate separation is decided by the sudden 
increase in the normal amplitude of the cantilever [20|. In our set-up the normal stiffness of the tuning fork- fibre 
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FIG. 8: The (a)Amplitude and (b)Phase of the tuning fork with separation at 12 KHz shear frequency. 
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FIG. 9; (a) The fraction change in drag force upon confinement from the bulk viscous drag, and (b)Relaxation time 
T with separation at 12 KHz shear frequency 



system is nearly infinite. We observe that the amplitude does not vary much after making a hard contact with the 
surface, this possibly means that the last layer of liquid on the surface is not drained out and provides a constant 
drag to the oscillating tip. The zero of the tip-substrate separation in our experiments is decided by considering the 
percentage change in the amplitude after the engagement. Using equation [6], the phase data is converted into real 
and imaginary part of damping coefficient. Equation(8) is used to plot relaxation time versus separation. Figure 9(a) 
shows a plot of relaxation time versus separation for a water at 12 KHz of shear frequency. The relaxation increases 
significantly compared to bulk water and is measurable under the shear at 12 KHz. It attains a value between 0.01 
and 0.1 ms when the film thickness is less than 5 nm, a value comparable to the one measured using the shear and 
normal measurements using AFM [19|, I20j . In these experiments the relaxation as well as drag force is measurable 
only below 1 nm of film thickness. This is possibly because the confining areas under AFM tip are much smaller 
compared to the present measurements and films with thickness above 1 nm offer immeasurably low change in drag 
force compared to the bulk viscous drag. Figure 9 (b) shows the relative change in the drag force when the tip is 
oscillating close to the substrate. It is instructive that the change in drag force offered by the nano-confined volume of 
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the liquid is only 2 to 8 percent of the bulk drag and will not have a discernible effect if the thermal noise in the force 
measurement is not low enough. By comparing the drag force and the relaxation time one can infer various regimes 
of non-Newtonian behaviour of the nano-confined liquid such as viscoelastic or shear thinning. The off resonance 
operation ensures that the measurements can be performed over a range of shear amplitudes and frequencies. 
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